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1. Introduction 

As the lightest quark not participating in the valence sector, the strange quark plays a special 
role in our understanding of the structure of the nucleon. In recent years, there has been a vigorous 
experimental effort to determine the strange contribution to the nucleon's elastic electroweak form 
factors. In particular, a number of experiments have sought to measure the strange electric and 
magnetic form factors via parity-violating electron scattering [|]], notably SAMPLE, A4, HAPPEX, 
and GO. Recent combined analyses [^, ^] find values for 0^(2^) and that are small and 

consistent with zero in the range of momenta so far explored. 

Also of interest is the strange axial form factor G\{Q^), to which electron scattering exper- 
iments are relatively insensitive. At present, the best constraints come from the two-decades old 
neutrino scattering data of the E734 experiment at Brookhaven [Q]. A recent analysis combin- 
ing these results with those of HAPPEX and GO, favors a negative value for G\{Q^) in the range 
0.45 < 2^ < 1.0 GeV^. In time, a more precise determination should be possible with results from 
the new generation of neutrino scattering experiments (e.g., BooNE, MINERvA). A special case is 
presented by the strange axial form factor at zero momentum transfer, G\ (0) = A^, which may be 
identified with the strange quark contribution to the spin of the nucleon. This quantity is of partic- 
ular importance, given the role sea quarks are thought to play in resolving the proton "spin crisis." 
In principle, it is accessible in deep inelastic scattering, where it is given by the first moment of 
the helicity-dependent structure function ^s{x). In practice, however, determining the first moment 
requires an extrapolation of the experimental data to small values of x, where uncertainties are less 
under control. There is some tension between the two most recent analyses from HERMES 0], 
which rely on different techniques; the former favors a negative value for As while the latter finds 
a result consistent with zero, within somewhat larger uncertainties. 

Unlike the strange electromagnetic and axial form factors, the strange scalar form factor 
G^^{Q^) is not directly accessible to experiment. At zero momentum transfer, this quantity corre- 
sponds to the strange scalar matrix element (A/^l^^lA/^). Often considered in relation to the nucleon- 
pion sigma term, it is an important parameter in models of the nucleon. It has also received atten- 
tion, however, for the role it plays in the interpretation of dark matter experiments. Many models 
of TeV-scale physics, including most supersymmetric extensions of the standard model, yield a 
dark matter candidate (e.g., neutralino) that scatters from nuclei via Higgs exchange. The Higgs 
predominantly couples to strange quarks in the nucleon, and so (A/^l^sjA/^) enters the cross-section 
through the quantity 

mMss\N) 

= ■ ^^-^^ 

As emphasized recently in [§, fis is poorly known at present and represents the leading theo- 
retical uncertainty in the interpretation of direct detection experiments. 

Strange form factors have been explored on the lattice in a number of past investigations, 
mainly within the quenched approximation. For a partial list of references, we direct the reader 
to [[Io|]. For some relevant interesting results presented at this conference, see [11, 12, pj| ]. 



2. Method 

We begin by defining the basic correlation functions from which we will extract the strange 
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form factors, starting with the usual two-point function for the proton. At momentum q, this is 
given by 

GP)(^?o;^) = (l + 74)"''Ie'^^"(/'^x,0^"(0,?o)). (2.1) 

.V 

Here = eabc{u^C'Y=,db)u'^ is the standard interpolating operator for the proton, with smeared 
quark fields, and (1 + 74) projects out the positive-parity state. Next, we define various three-point 
functions Gx {t,t' ,to;q^), where X = S,E,A correspond to the disconnected scalar, electric, and 
axial form factors, respectively. These are given by 

dpit/Mq) = {l + rAr^'£e"l-^^{pP{x,t)[ss{J^,t') - {ss{^ ,t'))]P''{d,to)) (2.2) 
for the scalar, 

G'i\tj',to-,q) = {\ + Y4r^l^e'''''{pP{x,t)[V4^,t')-{^^^^ (2.3) 
for the electric, and 

1=1. v,x' 

for the axial, where and are the point-split vector and axial currents. Note that we always 
employ the vacuum-subtracted value of the current, [J{x,t) — {J{x,t))], even though this is only 
strictly necessary when J is the scalar density, since the expectation value of the others vanish. 
Given finite statistics, however, and an inexact estimate of the trace, it is possible that using the 
vacuum-subtracted value gives reduced statistical errors. Empirically, we find that for the strange 
scalar and axial form factors at q^ = 0, the two approaches give indistinguishable results. At larger 
momenta, however, uncertainties for the vacuum-subtracted quantities are noticeably smaller. 
Inserting a complete set of states in Eq. ( p] ) yields the spectral decomposition 

G^'\t,to;q) = l.cMe-^"^'^'^'-"^ , (2.5) 

n 

where the coefficients are given by 

c„(^)=2(l + g)z2(^). (2.6) 

Here Z„(^) is the amplitude for annihilating the nth state with our interpolating operator (n = I 
corresponds to the nucleon itself). The momentum dependence in Z„ arises because we utilize 
smeared quark sources for the nucleon. Performing a similar decomposition for a generic three- 
point function, we find 

G^^\t,t',to;q) = '£jnm{q)e-"'"^'-''^e-^"'r^^(''-'»^ , (2.7) 

m,n 

where the coefficients j„„, may in general be expressed in terms of some combination of form 
factors. In particular, the correlation functions given in Eqs. QJl-lA) have been defined such that 
a single form factor enters the coefficient 71 1 for each case, according to 



iu{q)=2[ \ + j-^^]Z,{Q)ZMG'x{q') (2.8) 
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for X = S,E,A, where G^(^^) is the corresponding strange form factor of the nucleon. In terms of 
the coefficients c„ extracted from the two-point function G^^\t,tQ;q), this becomes 




;„(^) = Gi(^^)^/^( 1 + ^ )ci(d)ci(^). (2.9) 



Our general strategy will be to fit the correlation functions Gx {t,t' ,to;q) to Eq. (2.7), taking into 
account both the ground state nucleon and a single excited state. We may then extract the nucleon 
form factors from j\ i with input from the two-point function. In principle, one could also obtain 
form factors of the excited state from 722, as well as transition form factors from 712 and 721. In 
practice, however, we expect these to absorb the contributions of still higher states and trust only 
the ground state form factors to be reliable. 

The motivation for this approach will be discussed further as we present our results in the 
following section. Before turning there, however, we conclude this section with a few remarks on 
the calculation of the disconnected insertion. This represents a major challenge for the lattice since 
it requires the trace of the full Dirac propagator over color, spin, and spatial indices. Our approach 
for estimating this trace is described in [|l0||. Briefly, it involves inverting against a set of sources 
where each is nonzero on only a small number of sites, with these sites as widely separated as 
possible. In the parlance of the usual stochastic source method, this is equivalent to using a single 
noisy source with "extreme dilution." It provides a very accurate estimate of the trace, and since the 
small amount of contamination arising from off-diagonal terms is not gauge-invariant, the method 
is unbiased. On our 24^ x 64 lattice, we have chosen a dilution pattern where the smallest spatial 
separation between sites is 6\/3«s- With 864 (times 12 for color/spin) such sources per lattice, we 
obtain the trace over each of four time-slices. 

In order to construct the nucleon correlators, we also calculate 64 propagators at the light quark 
mass per lattice, with one originating from each time-slice. We utilize gaussian smearing for the 
Ught quark propagators. 



3. Preliminary results 

The results presented here were calculated on a 24^ x 64 anisotropic lattice, using an ensem- 
ble of 863 gauge configurations provided by the Hadron Spectrum Collaboration. These were 
generated with two dynamical flavors of unimproved Wilson fermions, with Mj^ 400 MeV and 
a,v = 0. 108(7) fm Sc?,. Given the anisotropy, this lattice has a relatively short extent in time, which 
has influenced our choice of method. Conventionally, one extracts the form factors by considering 
various ratios of the three- and two-point functions defined above. For example, at zero momentum 
transfer, one has 

R,{t,t'Mq'=0) = ^^ir^^GUq'), (3.1) 
G(2)(f,fo;0) 

for large time separations. As discussed in the previous section, we have chosen instead to fit 
the three-point function that appears in the numerator of this ratio directly. The reasons are two- 
fold. First, this allows us to avoid contamination from backward-propagating states, which are 
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problematic due to the short temporal extent of our lattice. At the same time, it allows us to 
explicitly take into account the contribution of (forward-propagating) excited states. 

Our fit of the scalar form factor at = is shown in Figure [I]. For the purpose of plotting, we 
have normalized our results by a fit to the two-point function. With this normalization, dominance 
of the ground state should manifest as a plateau at large times. The statistical error in the extracted 
form factor, GJ(0)iat = 2.70(38) is determined by jackknife over the full fitting procedure. We 
note that in this preliminary analysis, we have treated the system symmetrically by requiring that 
the separation between the source and current insertion equal that between the insertion and the 
sink, i.e., {t — t') = {t' — to). It should be possible to reduce the statistical errors further by lifting 
this constraint. We expect that this will prove especially advantageous for the results at nonzero 
momentum transfer presented below, since these require a fit to a sum of four, rather than three, 
exponentials. 

We have not yet attempted to quantify the systematic uncertainties. In addition to the usual 
systematics associated with the finite lattice spacing, finite volume, and the unphysically heavy 
mass of the light quarks, our approach is also sensitive to the choice of the fitting window in the 

(3) 

fit of the two-point function. Recall that in order to extract the form factors from G^^ {t,t' ,to;q), 
we must first determine the coefficients Cn{q) and masses/energies En{q) from a fit to G^^\t,tQ;q). 
Since we have access to a total of 863 x 64 = 55, 232 nucleon correlators, these tend to be very well- 
determined. The coefficients c„ are somewhat sensitive to the choice of fitting window, however, 
and since they multiply the form factor in Eq. (2^), this translates into a direct systematic error on 
the form factor, at perhaps the ten percent level. 

Given our determination of the strange scalar matrix element, for the renormalization-invariant 
quantity /ts, we estimate 

/n. = ^^^^^ = 0.51(8)(3), (3.2) 

where we have inserted the physical nucleon mass for M^. The second error is the uncertainty 
in relating this mass to the lattice scale, the first is statistical, and no other systematics have been 
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Figure 2: Results for the strange axial form factor Figure 3: Strange axial form factor as function of 
at = 0. momentum. 





Figure 4: Strange electric form factor as function Figure 5: Strange scalar form factor as function of 
of momentum. momentum. 



taken into account. The bare quark mass lUg is the naive mass appearing in the Wilson action, minus 
the critical mass, as determined from a (partially quenched) chiral extrapolation of the mass of the 
pseudoscalar. It is worth noting that the matrix element appearing in Eq. ( [3.2[ ) was calculated for a 
world with a 400 MeV pion. This is a consistent determination of the physical value only insofar 
as the matrix element (A/^j.wlA'^) is insensitive to the mass of the light quarks. As an alternative, we 
may work consistently in such a world by inserting our calculated value of the nucleon mass. In 
this case, the scale dependence drops out, and we find fjs = 0.44(7). 

Results for the strange axial form factor are shown in Figure]^ As for G^(0), we note that our 
result Ai' = G^(0) = —0.064(25) has not been renormalized and so may not be compared to the 
experimental results. It is noteworthy, however, that we find a value that is negative and distinct 
from zero at the level of 2.5a. In Figure^ we show the momentum dependence of G\{q^). 

Finally, in Figure |^ we present our preliminary results for the strange electric form factor as a 
function of momentum. Note that since the strange quark does not contribute to the electric charge 
of the nucleon, G|. (0) must vanish. For completeness, the momentum dependence of the strange 
scalar form factor is shown in Figure ||. 
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